Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



X-ray selected AGNs: The Clustering Properties 



E.Ghaffarsedeh 1 , J.Quenby 

1 Astrophysics Group, Blackett Laboratory, Imperial College, Prince Consort Rd, London SW7 2BZ, UK 



1 February 2008 

'66 : 

0\ • ABSTRACT 

G\ ■ 

Clustering properties of two samples of X-ray selected AGNs are compared with 
■ IR/Optical galactic samples. These are the Einstein Observatory Extended Medium 

Sensitivity Survey: EMSS (Maccacaro it et al. 1991), with flux limit of > 3.25(lCr 14 ) 
erg cm -2 sec -1 (0.5-2.0 keV), redshifts up to 2.83 and a sky coverage of 778 deg 2 ; 
and an all sky sample of 654 ROSAT detected radio quasars (Brinkmann et al. 1997), 
with a uniform flux limit of a few (10 -13 ) erg cm -2 sec -1 (0.1-2.4 keV) and redshifts 
up to 3.886. The cosmological dipole moment is calculated for these samples, yielding 
the values of bQ~ 0,6 (where b is the bias parameter, and Q Q is the present value of 
the cosmological density parameter). We find this factor to be 0.68 ± 0.22 for the 
qq ■ EMSS sample, and 1.02 ± 0.014 for the quasars sample. Comparing these values with 

' that of the IRAS 1.2 Jy sample (1.95 ± 0.096), suggests that these AGNs are less 

{SJ , clustered than the IRAS galaxies. For each sample there is a depth at which the 

t— I ' amplitude of the dipole growth curve saturates. This convergence radius (R CO nv) is 

t-H ; found to be - 300 - 400 h" 1 Mpc for the EMSS sample and - 450 - 500 h" 1 Mpc 

for the ROSAT detected quasars sample, which suggests that X-ray AGNs have a 
much deeper contribution to the Local Group (LG) motion compared to the galaxy 
^ (— | and cluster samples. Also the maximum amplitudes at saturation are found to be 

O 405 ± 103 km/sec for the EMSS AGNs and 617 ± 8 km/sec for the quasar sample, 

which suggests that the lower the flux limit of the sample (ie. the more sensitive the 
sample), the lower the maximum amplitude at the region of saturation, therefore the 
less clustered the sample. Using the properties of the angular and spatial correlation 
functions for these samples, we find: bi R As/b qU asars = 3.61 ± 2.7; bi R As/bEMSS = 
8.46 ±6.4 and b quasars /bEMSS = 2.35 ±0.2. Using Chi-squared minimisation, we have 
fitted the angular correlation functions using the standard power law model of the 
form: W(0) = A8~ s , and the spatial correlation functions using the power law model 
$2j ■ of the form: £(r) = (r/ro) -7 . For the angular correlation functions, we find A ~ 0.6 

' for the AGN auto-correlations and A ~ 0.2 — 0.3 for the AGN-IRAS galaxy cross- 

correlations. Also 5 is found to be 0.8 — 1.1 . In the spatial correlation analysis, we 
find r a ~ 3.5 — 5.1 h^Mpc and 7 ~ 1.4 — 1.7 . The bias factor resulting from the 
IR/Optical galaxies may not be the true representative of this parameter, due to the 
local nature of those samples, whereas AGN samples are deep/far enough to give a 
more accurate value of this parameter. 

Key words: AGNs - Clustering - Large scale structure of the universe - Cosmology: 
Observations 



o 



1 INTRODUCTION 

A way of investigating the large scale clustering properties of 
cxtragalactic objects, is to analyse their cosmological dipole 
moment. The anisotropy of the Cosmic Microwave Back- 
ground (CMB) radiation and the solar motion relative to 
the Local Group of galaxies (LG) implies that the LG has a 
peculiar velocity of ~ 600 km/sec towards (l,b) = (277° ,30°) 
(cf. Rowan- Robinson, Lynden-Bell, 1988). According to the 



linear perturbation theory (Peebles, 1980), the peculiar mo- 
tion of the LG is in the same direction as the gravitational 
acceleration arising from the distribution of mass outside of 
the LG. Attempts to trace this mass distribution respon- 
sible for the gravitational acceleration associated with the 
peculiar velocity of the LG have been made using the IRAS 
galaxies (Yahil, Walker & Rowan- Robinson, 1986; etc.), Op- 
tical galaxies (Lahav, 1987; etc.), X-ray clusters (Lahav et 
al., 1989) and X-ray AGNs (Miyaji & Bolt, 1990). 
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The spatial extend of the distribution of mass inhomo- 
geneities that cause the LG motion has been a matter of 
debate. The most probable cause for this motion as well as 
for the observed peculiar motions of other galaxies and clus- 
ters is gravitational instability. This is supported by the fact 
that the gravitational dipole (acceleration) of many differ- 
ent samples of extragalactic mass tracers is well aligned with 
the general direction of the CMB dipole (Rowan-Robinson 
et al., 1990; Strauss et al. 1992; etc.). 

However what still seems to be under discussion is the depth 
from which density fluctuations contribute to the gravita- 
tional field that shapes the LG motion. The largest such 
depth, is defined by the dipole convergence depth, R con v, 
which is that depth where the true gravitational acceleration 
converges to its final value. The apparent value of R CO nv dif- 
fers from sample to sample, normally in the range 40 to 100 
h" 1 Mpc (for galaxy samples), with a strong dependence on 
the sample's characteristic depth. Cluster samples, such as 
Optical Abell/ACO cluster sample is volume limited out to 
large enough depth (~ 240/i _1 Mpc) to allow a more reliable 
determination of R CO nv which was found to be ~ 160/i _1 
Mpc (Scaramella et al., 1991; Plionis & Valdernini, 1991; 
Branchini & Plionis, 1996). Therefore studies using clusters 
of galaxies as a probe, suggests that a significant part of 
the LG's motion is caused by structures at much greater 
distances than 100 h" 1 Mpc. Recently, this result has been 
suported using X-ray cluster samples, which are free of the 
various systematic effects from which the optical catalogues 
suffer (Plionis & Kolokotronis, 1998). 

In the following sections, we will be discussing the dipole 
moment for the EMSS AGNs and the ROSAT detected 
quasars sample, from which we will find the values of the 
bias parameters. Furthermore, to investigate the clustering 
of objects at small scales, we study their angular and spatial 
correlation properties. 

One of the most common statistical characterizations of the 
clustering properties of galaxies and other extragalactic ob- 
jects is the three-dimensional two point correlation function 
(Peebles, 1980), £(r), in which one uses information on red- 
shift, longitude and latitude for each object. Note that this 
is the Fourier transform of the power spectrum of the galaxy 
distribution (Peacock & Nicholson, 1991), where the power 
spectrum is the square of the absolute value of the Fourier 
transform of the density field of galaxies. An easier way to 
investigate these clustering properties, however, is the two- 
dimensional two point angular correlation function, w{6), in 
which we only need information on the longitude and lati- 
tude for each object. 

Later on, we will be discussing the results of the auto and 
cross-correlation functions between various samples, from 
which we can deduce the values of the relative bias param- 
eters. 



2 DESCRIPTION OF DATA 

The EMSS is a statistically complete (96%) and well-defined 
sample of 835 serendipitous X-ray sources (of which more 
than 400, are AGNs) detected in IPC images of the high 
Galactic latitude sky obtained with the IPC on board 
the Einstain Observatory. It has limiting sensitivities from 
5(10" 14 ) to 3(10" 12 ) erg cm" 2 sec" 1 in the 0.3-3.5 keV en- 



ergy band (typically > 3.25 (10~ 14 ) erg cm -2 sec" 1 in the 
0.5-2.0 keV energy band), redshifts up to 2.83 and a sky cov- 
erage of 778 deg 2 (for further details, see Stocke et al.,1991; 
Gioia et al.,1990; Maccacaro et al., 1991). 

The sample of 654 ROSAT detected radio quasars 
(Brinkmann et al., 1997), is an all sky survey with a uniform 
flux limit of a few (10~ 13 ) erg cm" 2 sec" 1 (0.1-2.4 keV) and 
redshifts up to 3.886. This sample has been compiled from 
the Veron-Cetty-Veron catalogue detected by ROSAT in the 
ROSAT All-sky survey (Voges, 1992), as targets of pointed 
observations, or as serendipitous sources from pointed obser- 
vations as available publicly from the ROSAT point source 
catalogue (ROSAT-SRC, Voges et al., 1995) 



3 A BRIEF REVIEW OF THE TESTS FOR 
COSMOLOGICAL EVOLUTION FOR XRAY 
SELECTED AGNS 

To test the evolution in cosmological objects, one normally 
investigates the logN-logS relation, the < V e /V a > test and 
the differential luminosity function. 

The logN-logS relation for the EMSS AGNs and the 654 
ROSAT detected quasars are shown in Figure 1. It appears 
that these AGNs have gone through a positive luminosity 
evolution with time, ie. they have either been more luminous 
or more common / numerous in the past. Also note that 
the flux at which the logN-logS curve stops increasing is an 
approximation for the flux limit of the sample. Errors have 
been estimated as pure Poissonian. Furthermore, using the 
maximum likelihood method, each curve has been fitted with 
a power law model . As can be seen, a steeper slope is found 
for the EMSS sample (-1.57 ± 0.1) compared to the 654 
quasar sample (—1.13 ± 0.36), from which it might appear 
that the EMSS AGNs have evolved more in time. However, 
considering the la errors, the slopes are consistent with -1.5 
which is the value found for non-evolving objects. 

Next, the differential X-ray Luminosity Function (XLF) 
for both the EMSS sample and the 654 quasar sample, are 
presented in Figure 2 and Figure 3. To obtain these, we have 
used the method of summing inverse accessed volumes, as 
discribed in Schmidt (1968); Maccacaro (1991) and Page et 
al. (1996). Therefore, the XLF is calculated as: 

i 

where, V s is the volume in the universe searched to find a 
particular object, and dL corresponds to the luminosity bins. 
Also the error bars are calculated using: 

JV 

(E^) 1/2 w 

Assuming only the luminosity evolution model (since this 
is more commonly accepted), a single power law model has 
been fitted to each individual XLF corresponding to different 
redshift shells and a different slope has been found for each, 
as shown on the plots. 

Note that it is more common to fit a LF with two power 
laws, but in this calculations, a single power law gave a rea- 
sonably good fit (The goodness of the fits can be tested using 
the 2D K-S test). Finding different slopes for each redshift 
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Figure 1. The differential logN-logS relation for the EMSS and 
the 654 ROSAT detected AGNs 
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logL[10 40 erg/sec : 0.5-2.0 ^eV] 

Figure 2. The differential XLF for the EMSS AGNs 



shell, suggests that these AGNs have evolved differently at 
different periods of time. 

The XLF is used to calculate the selection function, which 
is used in dipole and correlation function analysis. In dipole 
analysis, the weight corresponding to each object is pro- 
portional to the inverse of the selection function. Also in 
the analysis of the spatial correlation function, the random 
sample is selected so as to have the same selection function 
(ie. the same redshift distribution) as the original sample. 
In most common cases of flux limited extragalactic object 
catalogues, one has to take into account the effects of the 
consequent undersampling of the density field especially at 
large distances, where the radial selection functions rapidly 
decline. Assuming that the unobserved galaxies are spatially 
correlated with those included in catalogue, the usual proce- 
dure to correct for the missing population is to weight each 
observed object at a distance r by a factor oc , which is 
the reciprocal of the portion of the LF that can not be sam- 
pled at that distance due to the flux limit of the catalogue. 
Therefore, the weight becomes: 



: 4> 1 (n)i 



where: 



<j>(r\ 



$(L)dL 



is the selection function. 

$(L) is the luminosity function, L ri 



-l(r-) 



(3) 



(4) 



= 4nr Sum is the 



minimum luminosity that an extragalactic object can have 
in order to be visible at a distance r which is determined by 
the flux limit of the sample. L max is the maximum luminos- 
ity of such an object. 

Note that the average density, n, is calculated as: 



$(L)dL 



(5) 



log L (10 4U erg/sec : 0.5-2.0 keV) 

Figure 3. The differential XLF for the 654 ROSAT detected 
quasars 



where L m i„ is the minimum luminosity of the object. 
Note that the selection function is inversly proportional to 
r, and this acts as a compensating effect for the fact that 
poor sampling occurs at large redshifts. 



4 THE COSMOLOGICAL DIPOLE MOMENT 

In this section, we give a summary on how to calculate the 
amplitude of the dipole growth curve as a function of the 
radial distance. Basically, the dipole moment is defined as: 



i=N 

D = ^ Win 



(6) 



where r; is the unit vector pointing at the position of each 
object and N is the total number of such objects within the 
distance considered. The weight (wi) is calculated as shown 
in equation 3. 

Using the linear perturbation theory, and the above def- 
inition of the dipole, the relation between the observed pe- 
culiar velocity of an observer, Vp, and that predicted by its 
gravitational acceleration, g, as long as the two vectors are 
well aligned is: 
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V P (r) 



H (3 v-^ 4>~\ri)n 



Ann t-^ 



Pg{r) 



(7) 



where: (3 = — g— , (b is the bias parameter). 
Note that V p (r) decreases with increasing r, therefore the 
more distant objects have less contribution to the acceler- 
ation of the LG, but do continue to add to the peculiar 
velocity. 

From the above equation, the three components of the ac- 
celaration generated by each object on the LG can be cal- 
culated (V x , V y ,V z ). The force can then be smoothed out to 
guarantee linearity, where the function Smooth is defined as: 

1 



Smooth(r ,r min ) = 



(8) 



max(r 2 .r, r z min ) 

ie. each component is multiplied by the function smooth. 
We then sum over shells to build up the three cartesian com- 
ponents of the cumulative dipole (acceleration) at radius r. 
Finally we calculate the amplitude of the cumulative dipole 
at the same radius, using: 



V a 



2 

cum- 



4- V 



+ V 

I V C2 



(9) 



where V c 



Vr 



cum — x j v cum — y •, v cum 



. z are the three components 
of the cumulative acceleration. 

Note that the moments of the objects are affected by the 
discretness effects or the shot noise errors, which increase 
with redshift bacause of the rapidly declining selection func- 
tion. These effects introduce a variance in the dipole ampli- 
tude and direction. To calculate these effects, we can use the 
method of Strauss et al. (1992), Hudson (1993) and Bran- 
chini & Plionis (1996), as follows: 

i=N 

a 2 z = J2^(zi/\n\ 3 ) 2 (10) 

i=l 

Therefore: 

&3D = \f{(j'i + ^ + a\) (11) 

and the corresponding error along the line of sight (ID error) 
is: 

asD 



aiD 



3V2 



(12) 



5 THE RESULTS OF THE DIPOLE MOMENT 
ANALYSIS 

The resulting dipole growth curves for the EMSS sample is 
shown in figure 4. 

As can be seen, there is no signal observed at small 
depths (ie. < 40h~ l Mpc), indicating insignificant contri- 
bution to the LG motion from the nearby objects. This is 
expected, since there are not many objects in the EMSS 
sample with distances < 40ft~ 1 Mpc. The amplitude of the 
dipole growth curve saturates at rather large depths of 
> 300/i~ 1 Mpc. However from distances of about 300/i -1 
Mpc, the amplitude of the dipole increases only by small 
amounts. Therefore R CO nv ~ 300 — 400ft -1 Mpc, which is the 
depth from which these AGNs contribute to the acceleration 
of the LG via their gravitational attraction. Note that there 
are relatively large error bars, particularly at larger depths, 
which put limits to the accuracy of the value of Rconv ■ 




Figure 4. The Dipole Growth Curve for the EMSS AGNs 



The selected 352 IRAS galaxies hove the so 
i _ radial Selection Function os the EMSS AGNs 



i pure Poissonion) 




Figure 5. Comparing the Dipole Growth curves for the EMSS 
AGNs and the IRAS galaxies 



From this analysis, the value of ftjMSS^o ' 6 is found to be 
0.68±0.22. This depends on the present value of the cosmo- 
logical density parameter (cf. £l = 1 for the Einstein-de Sit- 
ter universe). For IRAS galaxies, R conv ~ 40 — 100ft -1 Mpc, 
therefore comparing this with the one for the EMSS sample, 
suggests a dependence between the maximum effective depth 
of the sample and the convergence radius, ie. the deeper the 
sample (in distances), the larger the value of Rconv 
To further compare the dipole moments of the EMSS AGNs 
with that of the IRAS 1.2 Jy galaxies, we have selected a 
sample of 352 IRAS galaxies with the same selection func- 
tion as the EMSS sample (ie. the same redshift distribution) 
and plotted the cumulative dipole growth curves for both 
samples, as shown in figure 5. 

In these analysis, we have plotted both curves with bias 
parameter of 1.0. As can be seen at almost all distances, 
the dipole moment of IRAS galaxies dominates that of the 
EMSS. It appears that IRAS galaxies have stronger contri- 
butions to the motion of the LG than the EMSS AGNs. 
In other words, IRAS galaxies are more clustered than the 
EMSS AGNs at scales of < 200/T 1 Mpc (this is compared to 
the underlying mass distribution), which is consistent with 
the values of the bias parameter (assuming the Einstein-de 
Sitter universe). 
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Figure 6. The Dipole Growth Curve for the 654 ROSAT detected 
radio quasars 



The selected 178 IRAS galaxies have the some red! 
Selection Function as the 654 ROSAT detected qua: 
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. IRAS 1 .2 Jy 
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Figure 7. Comparing the Dipole Growth curves for the 654 
ROSAT detected quasars and the IRAS galaxies 



Similarly, the dipole growth curve for the 654 ROSAT de- 
tected radio quasars, is shown in Figure 6. 
For this sample, R CO nv ~ 450 — 500ft~ 1 Mpc, which suggests 
that these AGNs contribute to the LG motion from much 
deeper distances than the IRAS or the EMSS samples. This 
can be explained as being due to the large redshifts of the 
objects in this sample. Also the amplitude at the saturation 
point, is found to be ~ 617±8fcm/sec, from which we obtain 



&654tj 



1.02 ± 0.01, which suggests that these 



AGNs are less clustered than the IRAS 1.2 Jy galaxies, for 
which b IRA stto 06 = 1-9 ±0.09. 

Finally, for a further investigations, we have plotted the 
dipole growth curves for this sample and a sample of 178 
IRAS 1.2Jy galaxies selected with the same selection func- 
tion, as shown in figure 7. 

Note that in these, we have assumed the bias parameter 
of 1.0 for both curves. This has created some overlapping 
between the two curves, suggesting a similar degree of clus- 
tering for both samples. 



6 THE ANGULAR AND SPATIAL 
CORRELATION FUNCTIONS 

The cross-correlation function £ a &( r ), between population 'a' 
and '6', is defined in terms of the joint probability 5P that 
a population 'a' object is in an infinitesimal volume element 
5V a and a population '&' in SVt separated from 8V a and SVt 
by r ah \ 



SP =< n a >< n b > [1 + £ab(rab)]SV a 6V b 



(13) 



where < n a > and < rib > are the mean number densities 
of population 'a' and respectively. 

If galaxies were distributed randomly in the universe, £(r) 
would be zero. If there is any distance scale over which ob- 
jects tend to cluster, £(r) > 0, and if there is a scale at 
which the objects avoid each other (anti-clustering), £(r) is 
negative within that scale. £(r) is calculated as: 

Nod (t) 



n R 
n D 



- 1 



(14) 



N DR (r) 

where Ndd and Ndr refer to the number of data-data and 
data-random pairs, respectively; no and ur refer to the 
mean densities of the real and random catalogues. 
Note that for the analysis of the spatial correlation function, 
the population of the random sample is chosen to be at least 
50 times more than the real sample. This random sample is 
selected so as to have the same selection function (ie. the 
same redshift distribution) as the real sample. 
This function has been calculated for galaxies (Davis & Pee- 
bles, 1983; Fisher et al., 1994; etc.), clusters of galaxies (Bah- 
call et al., 1984) and QSOs ( Shanks et al., 1987). 
Note that the errors are calculated as the weighted Poisso- 
nian: 



N, 



(15) 



N DR N dr 

A common feature of the spatial correlation functions, is 
that they can be approximated by a power law model (on 
small scales, r < 10 h _1 Mpc, where the galaxy distribution 
is characterized by strong nonlinear clustering) of the form: 

£(r) = (r/r )-< (16) 

where r a and 7 are the correlation length and the power 
law index respectively. For galaxies, r is of the order 5.4 
h _1 Mpc (Davis & Peebles, 1983) with 7 ~ 1.7 - 1.8 . This 
correlation length has been found to increase as the scale of 
the system increases (ie. a larger value of r is found for rich 
clusters and even a larger value is found for superclusters) . 
Bahcall & Soneira (1983) found r a ~ 25 h" 1 Mpc and 7 ~ 
1.7 for Abell clusters. 

To find the best power law fit, we have used the method 
of x 2 minimisation, and the errors on each parameter cor- 
respond to 8\ 2 = 1. Note that this method is valid only if 
the errors of the correlation function at a fixed separation 
are Gaussian distributed and if the values at different sep- 
arations are uncorrelated. The power law model is usually 
fitted to the data for r < 20 h _1 Mpc and since these separa- 
tions are typically much smaller than the size of the sample, 
the central limit theorem ensures that the distribution of 
errors will be Gaussian. 

According to the linear biasing theory (Kaiser, 1984), the 
two point correlation function between populations 'a' and 
'b' can be approximated by: 
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£,ab{r) = b a b b £ m {r) 



(17) 



where b a and &(, are the bias parameters (ie. the contrast 
enhancement factor of the distribution of tracers compared 
to the underlying mass distribution) for populations 'a and 
'b' respectively and £m(r) is the correlation function of the 
underlying mass density, 
b is defined as: 



\ p Jtracer 
(HE) 

\ p Jmass 



(18) 



Calculating the correlation function of various tracers indi- 
cates the relative bias parameters of these tracers (Lahav, 
Nemiroff et al., 1990). Also, comparing a cross-correlation 
function (£/x(r)) with an auto-correlation function (£ir(r))i 
gives extra relative biasing between two types of mass trac- 
ers: 

Zn(r) ~ bj (Ly > 

Having described the calculation of the spatial correlation 
function, let us now have a brief review of the angular cor- 
relation function. The angular correlation function can be 
expressed in terms of probability in a similar way as the 
spatial correlation function. If one considers two differen- 
tial elements of solid angle on the sky, dfli and (IQ2, then 
the joint probality, dP, that galaxies will occupy the two 
elements with an angular separation of 9 can be written as: 



dP = n 2 (l + w(6»))dn 1 dn 2 



(20) 



where n is the mean surface density of galaxies. Note that 
a random (Poisson) distribution of galaxies yields w(6) = 
for all 6. Also: 



N RR 



Ndd , 

w = — 1 



N DR 



(21) 



(22) 



Again, the errors are calculated as Poissonian errors. Finally, 
w(6) can be approximated by a power law of the form: 

-5 



w(o) = Ae~ 



(23) 



As before the best fit, can be found using the methods of \ 
minimisation, and the errors on each parameter correspond 
to S X 2 =l 



7 THE RESULTS OF THE ANGULAR AND 
SPATIAL CORRELATION FUNCTION 
CALCULATIONS 

The results of the various angular and spatial auto and cross 
correlation functions are shown in figures 8-15. 

The best power law fit to all the data analysed in this 
work are shown in solid lines. For further comparisons, we 
have also plotted power law fits from previous work. For an- 
gular correlation functions, these are: the fit found by Woods 
et al. (1997), studying the angular auto-correlation function 
for faint optical galaxies in high galactic latitude fields: w{9) 
= 0.508 6>~ ' 8 (shown as the dotted line); and the fit by 
Lidman et al. (1996), studying the two point angular auto- 
correlation function for a sample of optical galaxies selected 




1:Full line(EMSS Auto): A=0.6±0.03,d = 1 ,0±0. 1 
Reduced Chi-squared = 24.59 



2:Dotted line(Woods et ol.,1997): A=0. 508,6 = 0.8 
3:Dashed line(Lidman et al.,1996): A=0. 1 26,6=0.85 



t 1 — -++ - 

< t t t , • 



Figure 8. The EMSS AGNs Angular Auto-Correlation Function 



1}Full Nne[EMSS -ut- ■ i : 7=T -iv ■ ' • = 4±0.1 ^ 
Reduced Chi-squdred = 50/60 

2) Dot-Dash linedRAS 1 .2Jy(Fisher.1 994):ro=3.75,7 =1,' 

3) Dotted line:Galoxies(Groth&Peebles,1977):ro=4.7,7=1 

4) Doshed line:Abell clueters(Bdhcdll&Soneird,1 983): 

f(r)=360(rhr 1 '" 




r(tT 1 Mpc) 

Figure 9. The EMSS AGNs Spatial Auto-Correlation Function 



from two large format CCD surveys: w{&) = 0.126 <T ' 85 
(shown as the dashed line). 

For the spatial correlation functions, these are: The power 
law fitted by Fisher et al. (1994), calculating the spa- 
tial auto-correlation function for the IRAS 1.2 Jy galaxies 
(shown as the dot-dash curve): £(r) = (r/3.76)~ 1,66 ; a gen- 
eral fit found for the galaxy auto-correlation (Groth & Pee- 
bles, 1977), shown by the dotted line: £(r) = (r/4.7) -1 ' 8 ; and 
finally the fit found for the Abell clusters auto-correlation 
(Bahcall & Soneira, 1983), presented as the dashed line: 
£(r) = 360(r/i)- 18 . 

It is important to note that there is a linear dependence 
between the scale of the system and the correlation length 
(ie. the scale of clustering). r is found to be ~ 5h~ 1 Mpc 
for galaxies and ~ 25/i _1 Mpc for Abell clusters. 
Furthermore, we have compared the average amplitudes for 
various samples for different angular and distance separa- 
tions, as shown in Tables 1-4. 

Note that from all these comparisons, it appears that AGN 
/ quasar samples are less clustered than the galaxy samples, 
and this is consistent with the results of the dipole analysis. 
Note that AGNs are rare distance objects, and this puts lim- 
its to the accuracy of the correlation signals at small scales. 
In the spatial correlation functions, there are relatively large 
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1 ; Full line: El ISS-QS0):A=0.6±0. 1 4,6=1 .1 ±0.1 4 
Reduced Chi — squared = 4.08 



2:Dotted line(Woods et al.. 1 997): A=0. 508,6=0.8 
3:Dashed line(Lidman et al.,1996): A=0. 1 26,6 = 0.85 
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Figure 10. The EMSS AGNs-654 QSO Angular Cross- Figure 13. The EMSS AGNs-IRAS 1.2 Jy galaxies Spatial Cross- 

Correlation Function Correlation Function 
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1:Full line(IRAS-QSO): A = 0.3±0.02.6=0.8±0.32 
Reduced Chi— squared = 3.56 



2:Dotted line(Woods et al.,1997): A=0. 508,6 = 0. 8 
3:Dashed line(Lidrnan et al.,1996): A=0. 1 26,6=0.85 
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Figure 11. The EMSS AGNs-654 QSO Spatial Cross-Correlation 
Function 



Figure 14. The IRAS 1.2 Jy sample-654 QSO Angular Cross- 
Correlation Function 
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Reduced Chi-squared=2.19 
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3:Dashed line(l_idmcn et al.,1996): A=0. 1 26,6 = 0.85 
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Figure 12. The EMSS AGNs-IRAS 1.2 Jy galaxies Angular 
Cross-Correlation Function 



Figure 15. The IRAS 1.2 Jy sample-654 QSO Spatial Cross- 
Correlation Function 
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error bars at small scales (r < 107i - Mpc) and the data is 
consistent with zero. 

Finally, the relative bias parameters between various sam- 
ples (at r = 12.5 h _1 Mpc) are found to be: 

blRAS /b%54quasar = 3.6 ± 2.7 (24) 

biRAs/bsMSS = 8.5 ± 6.4 (25) 

be>54quasars/bEMSS = 2.4 ± 0.2 (26) 

which again suggests that AGNs / quasars are less clustered 
than the IRAS galaxies. 

8 DISCUSSION 

Having obtained different bias parameters for different sam- 
ples of AGNs, using the dipole moment and the standard 
methods of spatial and angular correlation functions, let us 
now briefly discuss the calculation of the spatial correlation 
function for a fractal. 

Analysis of galaxy and cluster correlations based on the 
concepts and methods of modern statistical Physics, led 
to the suggestion that galaxy correlations are fractals and 
not homogeneous to the limits of the available catalogs 
(Pietronero, 1987; Coleman et al., 1992). 
According to these analysis, galaxy structures are highly ir- 
regular and self-similar: all the available data are consistent 
with each other and show fractal correlations (with dimen- 
sions, D ~ 2) up to the deepest scales probed so far (eg. 
1000 h^ 1 Mpc). If galaxy distribution becomes really homo- 
geneous at a scale A within the sample in question (ie. if 
A <!C R e f f , where R e f / is the effective depth of the sample) , 
then one has r = A 2 1 ^ D ~ 3 \ For a fractal of dimension 
D ~ 2, therefore: r ~ A /2. If on the other hand, the fractal 
correlations extend up to the sample limits, then the result- 
ing value of r a has nothing to do with the real properties 
of the galaxy distribution, but it is fixed by the size of the 
sample. 

The expression of £(r) in the case of fractal distribution is: 

£M = ^p(-^-r 7 -i (27) 

where R s is the depth of the spherical volume where one 
computes the average density. 

The so-called correlation length, r , is a linear function of 
the sample size, R s , 

To = (^) 1/7 i? s (28) 

and hence it is a spurious quantity without physical mean- 
ing, but it is simply related to the sample size. Note that 
according to the above equation, there is a linear depen- 
dence between the scale/ size of the system and the corre- 
lation length , and this is consistent with the results from 
the standard power law model of the correlation function. In 
other words, as the size of the sample increases (ie.R s ), the 
value of the corrrelation length increases (cf. r ~ 5h~ 1 Mpc 
for galaxies, and r ~ 25h~ Mpc for Abell clusters). 
Finally note that if the distribution is fractal, even the dipole 
moment is affected, and its value depends on even very dis- 
tant objects, although weighted by 1/r 2 . 



9 SUMMARY OF CONCLUSIONS 

From the above analysis (the dipole moment and the corre- 
lation function), it appears that the EMSS AGNs are less 
clustered than the sample of 654 ROSAT detected quasars 
which itself is less clustered than the IRAS galaxies. This 
was numerically shown by the values of the bias parame- 
ters. 

From the dipole analysis, we deduced the values of R CO nv ~ 
300 - WOh^Mpc for the EMSS AGNs and ~ 450 - 
500h~ 1 Mpc for the quasar sample, from which it appears 
that these AGNs / quasars have a gravitational effect on the 
LG motion from much deeper distances than IRAS galaxies 
(Rconv ~ 40 - WOh" 1 Mpc) and Abell clusters (R conv ~ 
160/t -1 Mpc). Also the amplitudes of the dipole growth 
curves at the saturation region are found to be 405.4±103.36 
km/sec for the EMSS AGNs and 617.3 ± 8.35fcm/sec for 
the quasar sample, from which we find: bEMSS^o ' 6 = 
0.68 ± 0.22 for the EMSS sample, and b 65iquasars Q.~ °' 6 = 
1.02 ± 0.014 for the quasar sample. Comparing the maxi- 
mum amplitudes at the saturation region, suggests that the 
more sensitive the sample (ie. the lower the flux limit), the 
less clustered it is compared to the underlying mass distri- 
bution. 

Assuming that bQ,~ 6 ~ 1, which is supported by deep AGN 
samples, we conclude Q a ~ 1 if b ~ 1. However to obtain 
fi ~ 0.3 (Krauss, 1998), requires b significantly less than 1. 
From the correlation function analysis, it appears that the 
AGN / quasar samples are less clustered than the galaxies, 
which is consistent with the results of the dipole calculations. 
The values of the relative bias parameters were found to be: 
biR,As/be>54-quasar = 3.61 ± 2.7, biRAs/bEMSS = 8.46 ± 6.4 
and b6Si-quasars/bEMSS = 2.35 ± 0.2. However the valid- 
ity of these values is subject to the critisism that they are 
sample depth dependent. 
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e 

i.5° < e < 5° 



< WemSS > < W EMSS -(i54QSO(Cross) > < W IR AS > 

0.05 ±0.07 0.074 ±0.09 0.14 ±0.07 



< W Woods > < W Lidman > 
0.21 0.05 



Table 1. Comparing the average amplitude of the X-ray AGNs Angular Correlation functions with those of IR/Optical galaxy Auto- 
correlation functions 



Table 2. Comparing the average amplitude of the AGNs-galaxics Angular Cross-Correlation functions with those of IR/Optical galaxy 
Auto-Correlation functions 
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e 

9 < 1.5° 
1.5° < 6 < 5° 



< W EMS S-IRAS(Cross) > < W IRAS _ 654: Q S o(Cross) > < Wi R AS > < W W oods > < W Lidman > 

0.14 ±0.12 0.23 ±0.13 0.25 ±0.11 0.79 0.23 

0.14 ±0.06 0.05 ±0.07 0.14 ±0.07 0.21 0.05 
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r(h~ 1 Mpc) < S.EMSS > < £,EMSS-654QSO(Cross) > < ZlRAS > < (.galaxies > < (,Abellcluster s > 

r < 10 0.86 ±1.49 no signals! 1.15 1.78 39.4 

10 < r < 20 0.81 ±0.6 1.84 ± 1.2 0.11 0.13 2.96 

20 < r < 50 7.4 ±0.64 6.9 ± 0.5 0.03 0.03 0.7 

Table 3. Comparing the average amplitude of the X-ray AGNs Spatial Correlation functions with those of IR/Optical galaxy Auto- 
Correlation functions 



r(h 1 Mpc) < £,EMSS-IRAS(Cross) > < £/ RAS -654Q SO(Cross) > <£.IRAS> < {.galaxies > < £.Abellclusters > 

r < 10 0.58 ±0.34 10.3 ± 3.2 3.86 7.25 161.1 

10<r< 20 0.29 ±0.15 2.15 ± 0.58 0.1 0.13 3.03 



Table 4. Comparing the average amplitude of the AGNs-galaxies Spatial Cross-Correlation functions with those of IR/Optical galaxy 
Auto-Correlation functions 
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